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The current cosmological constraints on a dark matter axion are reviewed. We 
describe the basic mechanisms by which axions are created in the early universe, 
both in the standard thermal scenario in which axion strings form and in inflation- 
ary models. In the thermal scenario, the dominant process for axion production 
is through the radiative decay of an axion string network, which implies a dark 
matter axion of mass ma ~ 100 /ieV with specified large uncertainties. An infla- 
tionary phase does not affect this string bound if the reheat temperature is high 
Treh^/a Or, Conversely, for T^eh^/a, if the Hubble parameter during inflation is 
large H\ > /a; in both cases, strings form and we return to the standard picture 
with a ma ~ 100 /ieV dark matter axion. Inflationary models with /a > H\ > Treh 
face strong CMBR constraints and require 'anthropic misalignment' fine-tuning 
in order to produce a dark matter axion; in this case, some inflation models are 
essentially incompatible with a detectable axion, while others can be engineered 
to allow a dark matter axion anywhere in a huge mass range below ma ^ ImeV. 
We endeavour to clarify the sometimes confusing and contradictory literature on 
axion cosmology. 



1 Introduction 

The axion has consistently remained one of the most popular dark matter 
candidates ever since the early 1980's. Unlike many other more exotic particles, 
the axian's existence depends only on a minimal extension of the standard 
modeU'B which also solves one of its key difficulties — the strong CP problem 
of QCD. In the standard thermal scenario, cosmic axiona^are created primarily 
through the radiative decay of a network of axion stringaju formed in the early 
universe at t « 10~^^s; this network rapidly annihilates when the axion mass 
'switches on' at t « 10~'*s leaving no trace except for a background of cosmic 
axions. If these axions are to be the dark matter of the universe, then their 
mass must be about ma ^ 100/ieV (with significant uncertainties). □ On the 
other hand, if an inflationary phase in the early universe eliminates all axions 
strings with a low reheat temperature, then dark matter axions jare created by 
more quiescent mechanismgj and they can have any mass belowca ^ ImeV. 
Astrophysical constraints on the axion masaOH suggest that maJ;10meV, so 
a viable parameter window exists for the dark matter axion; the only doubt 
that remains is what it actually should weigh. ^ |_. 

Present large-scale axion search experimentaaM are focussing on a mass 
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range, ^ 1-10 /ieV, for a variety of historical and technological reasons. 
However, a number of proposed experiments may have the sensitivity_to look 
for the heavier dark matter axion of the standard thermal scenario,Ej which 
appears to have a better motivated mass prediction. This axion cosmology 
overview, then, keeps this missing matter quest to the forefront. Our primary 
focus is not on the viable axion mass range, but rather on cosmological pre- 
dictions of the mass of a dark matter axion, that is, an axion for which is 
of order unity. 

2 Axion properties and constraints 

The standard model of particle physics, based around the Weinberg-Salam 
model for electroweak interactions and QCD for strong interactions, has one 
significant flaw — the strong CP problem: non-perturbative effects of instantons 
add an extra term to the perturbative Lagrangian and the coefficient of this 
term, denoted 9, governs the level of CP violation in QCD. The absence of 
any such violation in all observed strong interactions imposes a constraint on 
9, the most stringent being 9 < 10"^" due to the absence of a neutron electric 
dipole moment £3 Since the value of 9 is effectively arbitrary, one is left with a 
severe fine tuning problem. p. 

The elegant solution of Peccei and Quinntl is to allow 9 to become a dy- 
namical field which relaxes toward the CP conserving value 9 = 0'| through the 
spontaneous breaking of a J7(l)-symmetry. Essentially, 9 becomes the phase 
of a complex scalar field = |(/!)| exp(i0), which takes a non-zero expectation 
value at the Peccei-Quinn scale /a, that is, it has a potential of the general 
form, 

= ^m' - fir + f.^ifj^ - cose) . (1) 

The azimuthal particle excitation in the ^-direction (around the bottom of 
the 'Mexican hat') is a pseudo-Goldstone bosono — the axion field a — fa9. 
Here, the second term in the potential gives a small mass ma to the axion, 
however, note that this mass ma(r) only 'switches on' through instanton effects 
at the QCD-scale Aqcd ~ 200MeV (that is, m^{T) « 0, T » Aqcd and 
ma{T) — nia, const., T « Aqcd). This ensures that the 0- field lies in the 
true CP-conserving minimum = at late times. Throughout this discussion, 
we consider the simplest = 1 axion which couples to a single quark and has 
a unique vacuum. Later, we shall comment on the A'^ > 1 axion with a cos NO 
potential, noting the implications of the additional vacuum degeneracy. 

"In fact, 9 7^ 0, since there are some CP violating weak interactions, and it has lp«pn 
shown that B > 10~^* in the standard model, although this is highly model dependent 
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The axion has an extremely smaU mass ma which is inversely proportional 
to the Peccei-Quinn scale /a, 



= 6.2 X 10~^eV 



( 



IQi^GeV 

/a 



) 



(2) 



The axion couplings to ordinary matter also scale as f^^^, specifically to nu- 
cleons ffaN = C^niT^/fa (N = n, p) and to photons gi,y = (a/27r/a)(C^ - 1.9) 
where Cn and C'-y are parameters of order unity. These expressions do not 
adequately encompass the model-dependencies of axion phenomenology which 
includes, for example, the DFSZ axion with a tree-level coupling to the electron 
and the KSVZ axion with none (refer to a review such as ref. [14]). 

Initially, it was supposed that the Peccei-Quinn scale was close to the 
electroweak phase transition /a ~ Iew; but an exhaustive search of accelerator 
data ruled out this possibility, implying /a^lO^GeV. However, the ensuing 
disappointment was short-lived because there is no phenomenological reason 
why the Peccei-Quinn scale /a could not be much higher, even up to grand- 
unification scales. Thus, the 'invisible' axion was born, an extremely light 
particle with almost undetectably weak couplings of order f~^. 

Accelerator limits on the new 'invisible' axion were soon superseded by 
astrophysical calculations of stellar cooling rates. Axions, being weakly cou- 
pled, can escape from the whole volume of a star and, in certain parameter 
ranges, this can exceed the usual heat loss mechanisms through ccpvection and 
other surface effects. The oft-quoted 'red giant limit' on the axionlj is actually 
slightly weaker than that obtained from old globular cluster stars|0 TOa^0.4eV 
or /a?i2 X 10^. However, more stringent constraints come from supernova 
1987a; the duration of the observed neutrino signal, associated with nascent 
neutron star cooling, depends on the axion- nuclcon coupling (jaN- Originally, 
a strong bound rria^ ImeV was proposed,Ll but this did not adequately reflect 
the model-dependence of axion phenomenology. More conservative estimates 
in realistic axion models yield the overall boundQ 



Note that significant uncertainties remain for all these astrophysical constraints 
and recent reviews should be consulted such as refs. [7,8]. Interestingly, as- 
trophysics seems to have no impact in any cosmological scenario on the viable 
mass range of a dark matter axion. 



TO; 



< lOmeV, 



/a > lO^GeV. 



(3) 
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3 Standard thermal axion cosmology 

The cosmology of the axion is determined by the two energy scales /a and 
Aqcd • The first important event is the Peccei-Quinn phase transition which is 
broken at a high temperature T ~ lO^GeV by (||). This creates tha axion, 
at this stage effectively massless, as well as a networJtxif axion stringstj which 
decays gradually into a background cosmic axions.Hu One can, in principle, 
engineer models in which an inflationary epoch interferes with the effects of 
the Peccei-Quinn phase transition, but we shall deal with this case in the next 
section, focussing first on this 'standard thermal scenario'. At a much lower 
temperature T Aqcd after axion and string formation, instantnn effects 
'switch on' the axions acquire a small mass, domain walls formlHI between 
the stxingstS and the complex hybrid network annihilates in about one Hubble 
time .113 In the cold dark matter scenario of interest here, the remaining axions 
redshift and eventually dominate the energy density of the universe at equal 
matter-radiation t^q. 

There are three possible mechanisms by which axions are produced in the 
'standard thermal scenario': (i) thermal production, (ii) axion string radiation 
and (iii) hybrid defect annihilation when T = Aqcd- Axions consistent with 
the bound (0) decouple from thermal equilibrium very early at a temperaturelij 
Td<;4 X lO^cV; their subsequent history and number density is analogous to 
the decoupled neutrino, except that unlike a lOOeV massive neutrino, ther- 
mal axions cannot hope to dominate the universe with ma^lOmeV. Thermal 
production can still continue down to about Aqcd but this contribution is 
subdominant in light of the astrophysical bound (H) . We now turn to the two 
dominant axion production mechanisms, but first we address an important 
historical digression. 

3.1 Misalignment misconceptions 

The original papers on axions suggested that axion production primarily oc- 
curred, not through the above mechanisms, but instead by 'misalignment' ef- 
fects at the QCD phase transition. eIlj As applied to the thermal scenario, 
this turned out to be a very considerable underestimate, but it is useful to 
recount the argument and its flaws because they are often repeated in the 
literature. Before the axion mass 'switches on', the axion field takes random 
values throughout space in the range to 2tt. However, afterwards the true 
minimum becomes 9 — 0, so the field in the 'misalignment' picture begins to 
coherently oscillate about this minimum; this homogeneous mode corresponds 
to the 'creation' of zero momentum axions. Given an initial rms value 9i for 
these oscillations, it is relatively straightforward to estimate the total energy 
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density in zero momentum axions and compare to thejaEeaent mass density of 
the universe (assuming a flat il = f FRW cosmology) :Q'E2rE^ 

O., . 0.2 A^- .ft,.,) (^)"^ 2A.-.?f,.,,(i5^)" " 

(4) 

where A « accounts for model-dependent axion uncertainties, as well as 
those due to the nature of the QCD phase transition, and h is tbe rescaled 
Hubble parameter. The function i{6) is an anharmonic correctioncJ for fields 
near the top of the potential close to unstable equilibrium 9 ~ that-is, with 
f(0) = 1 at the base « and slowly diverging approximately asESl i[9\.^ 
— ln(|7r — 0|) for — > tt. If valid, the estimate would imply a constraintBEJ 
/aJilO^^GeV, TOa^5/xeV for the anticipated thermal initial conditions with 
6*1 (suggested to be 6; = tt/VS). 

In the thermal scenario, however, the expression (^ has three serious short- 
comings: (i) The axions are not 'created' by the mass 'switch on' at i = tgcD, 
they are already there with a specific momentum spectrum g{k). Dynami- 
cal mechanisms prior to this time — thermal effects, string radiation and field 
realignment — will have created this axion spectrum. Because correlations can- 
not be established acausally, there will be a lower cutoff to the axion momenta 
set by the horizon, fc"?^^ ^ t; only inflationary models can have truly zero mo- 
mentum axions. The actual axion number obtained from g(k) is clearly much 
larger than an rms average in (|) which ignores the true particle content, (ii) 
Secondly, this estimate was derived before much stronger topological effects 
were realized, notably the presence of axion strings and domain walls. These 
nonlinear effects complicate the oscillatory behaviour of 9 considerably, imply- 
ing that the homogeneous estimate (^) is poorly motivated, (iii) Finally, from 
a more general perspective, the expression (^ with /a » lO^^GeV seems to 
imply that axion mass 'switch on' changes the topology of the universe, tak- 
ing us to an r2 > 1 closed model! Of course, this is not what happens, since 
the extrapolation breaks down and in reality the axions merely dominate the 
universe much earlier than tcq, we end^jip with an fia — *■ 1 model with an 
unacceptably low baryon-to-axion ratio£3 

3.2 Axion string network decay 

Axions and axion strings are inextricably intertwined. Like ordinary supercon- 
ductors or superfluid ^He, axion models contain a broken [/(l)-symmetry and 
so there exist vortex-line solutions. Combine this fact with the Peccei-Quinn 
phase transition, which means the field is uncorrelated beyond the horizon, 
and a random network of axion strings must inevitably form. An axion string 
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corresponds to a non-trivial winding from to 2tt of the axion field 9 around 
the bottom of the 'Mexican hat' potential (|^). It is a global string with long- 
range fields, so its energy per unit length fi has a logarithmic divergence which 
is cut-off by the string curvature radius R^t, that is, 

M«2^/fln(t/5), (5) 

where the string core width is (5 /^^^ (taking A « 0{1) in (|l|)). By observing 
the logarithmic dependence in (j^), the axion string might be interpreted as a 
delocalised 'cloud of energy', but on cosmological scales it is anything but non- 
local; if we have a string stretching across the horizon at the QCD temperature, 
then \n(t/d) 65 and over 95% of its energy lies within a tight cylinder 
enclosing only 0.1% of the horizon volume. To first order, then, the string 
behaves like a local cosmic string, a fact that can be established by a precise 
analytic derivation and careful comparison with numerical simulations. || 

After formation and a short period of damped evolution, the axion string 
network will evolve towards a scale-invariant regime with a fixed number of 
strings crossing each horizon volume (for a cosmic string review see ref. [33]). 
This gradual demise of the network is achieved by the production of small loops 
which oscillate relativistically and radiate primarily into axions. The overall 
density of strings splits neatly into two distinct parts, long strings with length 
£ > t and a population of small loops £ < t, 



P = Poo + PL- (6) 

High resolution numerical simulationsBS confirm this picture of string evo- 
lution and suggest that the long string density during the radiation era is 
Poo ~ 13/i/t^. To date, analytic descriptions of the loop distribution have used 
the well-known string 'one scale' model, which predicts a number density of 



''The fact that a global string behaves like a local cosmic string, with some additional 
radiative damping, can be mathematically proven by renormalizing tlw-|Stiifi|igncptprgy and 



showing that it obeys exactly the same Nambu equations of motion.Ej'BjoEj'Cfl Numer- 
ical simulatirttift-nan then be used to establish the quantitative veracity of this analytic 
description.CicjEj However, great care must be taken in extrapolating string field theory 
simulations to cosmological scales, because the numerical results inevitably have a limited 
dynamic range which is over 10^® times smaller! For example, on these small numerical 
length scales global strings are in a regime where radiative damping is more than an order of 
magnitude stronger. If such numerical results are simply extrapolated to cosmology without 
a precise quantitative mathematical basis for describing the string motionj— tJtefi it might 
appear that axion strings radiate rapidly with a high frequency spectrunj: 
sensus in the literature, however, weighs heavily against such a conclusion £2lt 
Generically, cosmic axion strings radiate like a classical source with a spectrum weighted 
towards the lowest available harmonics. 
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loops defined as ij£n{£, t) i 
by! 



Pl{(-i t)di in the interval £ to £ + d£ to be given 
4ai/2(i + ^/^)3/2 



{£ + K<)5/2i3/2 



(7) 



where a is the typical loop creation size relative to the horizon and k « 
65/[27rln(<:/(5)] is the loop radiation rate. Once formed at t = to with length £o, 
a typical loop shrinks linearly as it decays into axions £ — £o ~ K(t ~ to) . The 
key uncertainty in this treatment is the loop creation size a, but compelling 
heuristic arguments place it near the radiative backreaction scale, a ^ k. 

String loops oscillate with a period T — £/2 and radiate into harmonics 
of this frequency (labelled by n), just like other classical sources. Unless a 
loop has a particularly degenerate trajectory, it will have a radiation spectrum 
Pn oc with a spectral index q > 4/3, that is, the spectrum is dominated by 
the lowest available modes. Given the loop density (0), we can then calculate 
the spectral number density of axions. 



!-(! + - 

K 



(8) 



which is essentially independent of the exact loop radiation spectrum provided 
q > 4/3. From this expression we can integrate over lu to find the total axion 
number at the time tqcDi that is, when the axion mass 'switches on' and the 
string network annihilates. Subsequently, the axion number will be conserved, 
so we can find the number-to-entropy ratio and project forward to the present 
day. Multiplying the present number density by the axion mass ma will give 
us the overall string contribution to the density of the universes 



llAh- 



1 + ^ 

K 



3/2 



/a 



10i2GeV 



1.18 



(9) 



For the expected ratio a/n « 1, the axion contribution arising from direct 
radiation from long strings is an order of magnitude smaller than f2a/, but it 
slightly reinforces (||) U We note also that fia/ is well over an order of magnitude 
larger than the 'misalignment' estimate (^). The key additional uncertainty 
from the string model is the ratio a/K, which should be clearly distinguished 
from particle physics and cosmological uncertainties inherent in A and h, and 
appearing in other estimates of fia- With a Hubble parameter near h = 0.5, 
the string estimate tends to favour a dark matter axion with a mass ma ~ 
lOO/xeV, as we shall discuss in the conclusion. 
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3.3 Axion mass 'switch on' and hybrid deject annihilation 



Near the QCD phase transition the axion acquires ajroass and network evolu- 
tion alters dramatically because domain walls form.ll30 Large field variations 
around the strings collapse into these domain walls, which subsequently begin 
to dominate over the string dynamics. This occurs when the wall surface ten- 
sion a becomes comparable to the string tension due to the typical curvature 
a ~ /i/t. The demise of the, hybrid string- wall network proceeds rapidly, as 
demonstrated numerically.Eai23 The strings frequently intersect and intercom- 
mute with the walls, effectively 'slicing up' the network into small oscillating 
walls bounded by string loops. Multiple self-intersections will reduce these 
pieces in size until the strings dominate the dynamics again and decay contin- 
ues through axion emission. 

An nrder-of-magnitude estimate of the demise of the string-domain wall 
networkt3 indicates that there is an additional contribution 

n,dw-0(l)A/.-^(^^^^j . (10) 

This 'domain wall' contribution is ultimately due to loops which are created 
at the time ^ iqcD- Although the resulting loop density will be similar to 
(^, there is not. the same accumulation from earlyptimes, so it is likely to 
be subdominantcl relative to (^). More recent workJSZl questions this picture 
by suggesting that the walls stretching between long strings dominate and 
will produce a contribution anywhere in the wide range fia.dw (l~44)r2a,£; 
however, this assertion definitely requires further quantitative study. Overall, 
the domain wall contribution will serve to further strengthen the string bound 
(^ on the axion. 

We note briefly that it is possible to weaken any axion mass bound through 
catastrophic entropy production between the QCD-scale and nucleosynthesis, 
that is, in the timescale range_lQr*s^tont ^ lO^^s. The idea has been sug- 
gested in a number of contexts,Lj'Ej but these usually involve the energy density 
of the universe becoming temporarily dominated by an exotic massive particle 
with a tuned decay timescale. 



4 Inflationary scenarios 



The relationship between inflation and dark matter axions is enigmatic. Its 
significance depends on the magnitude of the Peccei-Quinn scale /a relative 
to two key inflationary parameters, (i) the reheat temperature of the universe 
Tj-eh at the end of inflation and (ii) the Hubble parameter Hi as the observed 
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universe first exits the horizon during inflation. Inflation is irrelevant to the 
axion if T^ch^fa. because, in this case, the PQ-symmetry is restored and the 
universe returns to the 'standard thermal scenario' in which axion strings form 
and the estimate pertains. Even if the reheat temperature is low Tich < fa, 
however, inflation will again be irrelevant if Hi > because axion strings will 
form towards the end of inflation and we will return to a modified 'standard 
thermal scenario' — as we shall discuss. 

Thus, for inflation to impact the viability of a dark matter axion we must 
have /a^-ffi ^Tl-eh, that is, inflation must occur at a low energy scale or we 
must have an extremely light axion. Inflation, in this case, essentially makes 
no prediction as to the axion mass. Like much in axion cosmology, these facts 
are not widely appreciated, so they deserve some case-by-case unravelling. 

^.1 Hi < /a.' Anthropic misalignment and quantum fluctuations 

In an inflationary model for which /a > Hi > Treh, the ^-parameter or ax- 
ion angle will be-set homogeneously over large inflationary domains before 
inflation finishes .c3 In this case, the whole observable universe emerges from a 
single Hubble volume in which this parameter has some fixed initial value 9i. 
Because the axion remains out of thermal equilibrium for large /a, subsequent 
evolution and reheating does not disturb 9{ until the axion mass 'switches on' 
at T ~ Aqcd- Afterwards, the field begins to oscillate coherently, because it is 
misaligned by the angle 9i from the true minimum 9 = Q. This homogeneous 
mode corresponds to a background of zero momentum axions and it is the 
one circumstance under which the misalignment formula (jj) actually gives an 
accurate estimate of the relative axion density f2a- 

By considering the dependence r2a,h oc 0p in (^), we see that infiation 
models have an intrinsic arbitrariness eiyen by the different random magnitudes 
of 9i in different inflationary domainsO While a large value of /a » lO^^GeV 
might have been thought to be observationally excluded, it can actually be 
accommodated in domains where 9^ « 1. This may seem highly unlikely but, 
if we consider an infinite inflationary manifold or a multiple universe scenario 
including 'all possible worlds', then life as we know it would be excluded from 
those dcpaains with large 9i = 0{1) because the baryon-to-axion ratio would be 
too low.c3 Thus, accepting this anthropic selection effect, we have to concede 
that axions could be the dark matter f2a ~ 1 if we live in a domain with a 
'tuned' 0-parameter: ^ 

9i « A-'/^h ( ^ « 0.3A-1/2/, f^!^y ' . (11) 

'^This is not quite in the spirit of the original motivation for the axion! 
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For 9i « 0{1), this suggests an axion with Wa ^ 5/icV {h = 0.5), though 
actually inflation makes no definite prediction fiMfi ( pi] ) beyond specifying 
"^a^ 10~^eV. But even this restriction is not valid;Lll if we observe (^ carefully 
we see that we can also obtain a dark matter axion for higher Wa by fine-tuning 
9i near tt. The anharmonic term i{9) with an apparent logarithmic divergence 
allows r^a ~ 1 foiS \tt - 9i\ - exp[-1.5(/a/10"GeV)i-2] with /a<10"GeV, 
that is, for a much heavier axion. 

This simple picture is altered considerably when we include quantum ef- 
fects. Like any minimally coupled massless field during inflation, the axion will 
have a spectrum of quantum excitations associated with the Gibbons-Hawking 
temperature T ~ H/2t:. This implies the field will acquire fluctuations about 
its mean value 9i of magnitude 

69 - i//2^/a , (12) 

giving an effective rms value 9^ff = {9i + 69)^. Even if our universe began in 
an inflationary domain with 9i = 0, there will be a minimum misalignment 
angle set by 69; this implies that we cannot always fine-tune 9i in (|ll[) such 



that rici<l. Taking 9i > 69, we can compare equations (11) and (|12D to 
constrain Hi roughly as H < f^{fJlO'^'^GeV)-° '^ for U ^ lO^eV (h = 0.5). 
For larger axion masses, we instead require [tt — > 59, so we can similarly 
constrain iJ </aexp[-1.5(/a/10"GeV)-i-2] for /a<10"GeV. We illustrate 
this restriction schematically in fig. 1, showing that dark matter axions can be 
compatible with suitably constructed inflation models for any ma^lmeV (see 
ref. [22]). 

The quantum fluctualipns 69, in turn, will result in isocurvature fluctua- 
tions in the axion densityEJ given by (jj), that is, Spa oc 9^^ w 9f + 29^59 + 59'^. 
Such density fluctuations will also create cosmic microwave background (CMB) 
anisotropics which are strongly constrained, 6T/T ~ 6p/p ^ 2^0/6'cff ^ 10^^. 
Combining the fluctuation ( |l2|) with the ^-requirement for a dark matter ax- 
ion ( |Tl| ) wp nhta.in another strong constraint on the Hubble parameter during 
inflationCJOO 



0.4 



Here, Hi is the Hubble parameter as fluctuations associated with the time 
when the microwave anisotropics first leave the horizon some 50-60 e-foldings 
before the end of inflation. It is related to the vacuum energy at this time 
by Vi^'^ ~ {Hi X 10i®GeV)i/2. We conclude from (|l|) that inflation in the 
/a > > Troh regime must have a small Hubble parameter Hi or the dark 
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Figure 1: Inflation and the dark matter axion. Constraints on the Hubble parameter 
Hi during inflation and the Peccei-Quinn scale /a and axion mass ma. Only the white areas 
are compatible with a dark matter axion 1 in simple inflation models. For H > /a, the 

string bound (b|) pertains and, for H < /a, consistency is due to 'anthropic flne-tuning'. 



matter axion must be extremely light, q as summarized in fig. 1. 

Many inflationary models are not very consistent with the requirement 
(|l3|). For example, the simplest chaotic inflation models have Hi ^ 10^^-10^'* 
GeV, implying a truly invisible axion with ma ^ lO^^^eV and /a beyond the 
Planck mass! However, inflation now has many guises and it is possible to engi- 
neer models with Hubble parameters anywhere in the range H ~ 10^-10^''GeV; 
here, the lower bound is set by supersymmetry considerations and the upper 
bound by limits on adiabatic density fluctuations (see, for example, ref. [45]). 
New inflation can haveea iJ<10^GeV with a detectable dark matter axion, 

''Notethat if we were not considering a specifically dark matter axion satisfying (pj|), the 
bound (|l3|) is considerably weaker. 
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but the low H accentuates its severe initial condition difficultiesEJ (see re- 
cent axion inflation reviews such as ref. [47]). Another way to circumvent the 
CMB constraint (|l^) is in more complicated_LiiEationary scenarios in which 
the Peccei-Quinn scalar fielcL^ is the inflatonocj or else interacts with it, as 
in a hybrid inflation modelo It is then possible to have \4>\ » /a as the in- 
flaton rolls down to the vacuum, so the quantum fluctuations in (|l^) will be 
suppressed H by the factor fs./\<t>\- These inflationary caveats and the lack of 
any definite axion mass prediction has led LindecJ to conclude that "the best 
resolution of the uncertainties would be given by an experimental measurement 
of the axion mass. ... It may also help us to choose between different versions 
of inflation and axion models." ! 



4-.2 Hi^fi^: Axion string creation during inflation 

Even for a low reheat temperature Treh < /a, one can envisage inflation models 
with a large Hubble parameter during inflation Hi ^ fa.. The quantum fluctu- 
ations (|lj) in this case are sufficient to take the Peccei-Quinn field 4> over the 
top of the potentialcj leaving large spatial variations and non-trivial windings 
in (j). We can interpret this as the Gibbons-Hawking temperature 'restoring' 
the PQ-symmetry Ton^/a- As inflation draws to a close and H falls below 
/a, these fluctuations will become negligible and a string network will form 
along lines where ^ = 0. Provided inflation does not continue beyond this 
point for more than about another 30 e- foldings, we will effectively return to 
the 'standard thermal' scenario in which axions are produced by a decaying 
string network. It is effectively as if Treh > /ai-^o such low reheat inflation 
models are compatible with a dark matter axionE3, ~ lOO/ieV. This leaves 
the additional small upper window for the axion shown in fig. 1. 

More problematic is the case where H ^ /a, or just below, for an ex- 
tended period, say more than n > 30 e-foldings. In this case, string formation 
may be rare, but the field will still be driven around the Mexican hat po- 
tential by quantum fluctuations ending up with an average total fluctuation 
50 ^ y^iJ/27r/a<ll. This will imprint non-trivial windings into the axion 
field which will appear as exponentially large domain wallscj when the axion 
mass 'switches on'. Depending on their sizCpand distribution these are po- 
tentially incompatible with CMB constraints. C3 A number of other complex 
hybrid scenarios can arise when Hi fn, as reviewed in ref. [48]. 

Up to this point we have only considered the simplest axion models with a 

'^Such mechanisms might be thought to increase the inflationary upper bound rria ^ ImeV 
on therdprk matter axion, but the efli'ect is not dramatic because of the exponential fine- 
tuningE3 
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unique vacuum iV = 1, so what happens when > 1? In this case, any strings 
present become attached to N domain walls at the QCD-scale. Such ajiintj 
work 'scales' rather than annihilates, and so it is cosmologically disastrous .ESej 
Consequently, only /a > -ffi > Troh inflationary scenarios are acceptable for 
N > 1 and, for these, constraints like (Q) can be suitably adapted by merely 
rescaling /a ^ f^/N. 

5 Conclusions 

We have endeavoured to provide an overview of axion cosmology focussing on 
the mass of a dark matter axion. First, the cosmological axion density has been 
calculated in the standard thermal scenario by considering the dominant contri- 
bution from axion strings. In this case there is, in principle, a well-defined cal- 
culational method to precisely predict the mass ma of a dark matter axion. For 
the currently favoured value of the Hubble parameter {Hq k, 60kms^^Mpc~^), 
the estimate predicts a dark matter axion of mass 

TOa ~ 200 /xeV , /a ~ 3 X 10i°GeV , (14) 

where the significant uncertainties from all sources approach an order of mag- 
nitude. The key uncertainty in this string calculation is the parameter ratio 
a/ft, that is, the ratio of the loop size to radiation backreaction scale. Here, like 
most authors we have assumed a/n w 1, however, reducing this uncertainty 
remains a key research goal, though a technically difficult one. 

Secondly, we have reviewed inflationary axion cosmology showing that 
(i) many inflation models return us to the standard thermal scenario with 
ma ^ lOO/icV, (ii) some inflation models are essentially incompatible with 
a dark matter axion and (iii), because of the possibility of 'anthropic fine- 
tuning', other inflation models can be constructed which incorporate a dark 
matter axion mass anywhere below ma^ ImeV. 

We conclude that, while a dark matter axion might possibly lurk anywhere 
in an enormous mass range below ma^lmeV, the best-motivated mass for 
future axion searches lies near ma ~ 100/xeV, a standard thermal scenario 
prediction which is also compatible with a broad class of inflationary models. 
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